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Abstract. Correlated basis function perturbation theory and the formalism of cluster expansions 
have been recently employed to obtain an effective interaction from a state-of-the-art nucleon nu- 
cleon potential model. The approach based on the effective interaction allows for a consistent de- 
scription of the nuclear matter ground state and nucleon-nucleon scattering in the nuclear medium. 
This paper reports the the results of numerical calculations of different properties of nuclear and 
neutron matter, including the equation of state and the shear viscosity and thermal conductivity 
transport coefficients, carried out using the effective interaction. 
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INTRODUCTION 

The theoretical description of neutron stars requires a quantitative understanding of both 
equilibrium and non equilibrium properties of cold nuclear matter at high density. While 
equilibrium properties, e.g. the equation of state (EOS), are generally obtained from 
realistic dynamical models, strongly constrained from nuclear systematics and nucleon- 
nucleon scattering data, the studies of non equilibrium behavior, based on the solution of 
the linearized Boltzmann equation, often resort to oversimplified models of the nucleon- 
nucleon (NN) interaction. 

It has been recently suggested iQlfllll] that the many body formalism based on corre- 
lated wave functions and cluster expansion techniques provides an ideal framework for 
the unified treatment of a variety of nuclear matter properties, including the EOS, the dy- 
namic response, the spin susceptibility and the transport coefficients. The main element 
of the approach of Ref. HI, Q, 1] is a well behaved effective interaction, derived from 
a state-of-the-art model of nuclear dynamics, suitable for use in standard perturbation 
theory. 

In this paper, after reviewing the underlying assumptions of nonrelativistic nuclear 
many body theory (NMBT), we briefly discuss the derivation of the effective interaction 
within the Correlated Basis Function (CBF) approach, and focus on its application 
to the description of NN scattering in the nuclear medium. The results of numerical 
calculations of the shear viscosity and thermal conductivity of pure neutron matter, 
carried out using the in-medium NN cross section within the framework of Landau 
theory of normal Fermi liquids, are also reported. 



FORMALISM 



The paradigm of Nuclear Many Body theory 

Within Nuclear Many Body Theory (NMBT), nuclear matter is viewed as a collection 
of pointlike protons and neutrons, whose dynamics are described by the hamiltonian 
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m and p, being the nucleon mass and the momentum of the i-th particle, respectively. 
The NN potential v ;/ reduces to the well known Yukawa one-pion exchange potential at 
large distances, while its short and intermediate range behavior is determined by fitting 
the available experimental data on the two-nucleon system (deuteron properties and ~ 
4000 NN scattering phase shitfs). The state-of-the art parametrization of Wiringa et al 
0], generally referred to as Argonne vis, is written in the form 
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where 

aff 6 = [\,{p v o j ),S ij \®[\,{Xr'C j )], (3) 
Oi and T, are Pauli matrices acting in spin and isospin space, respectively, and 
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Sij = — (Oi ■ rtj) (a j ■ Tij) -(o r Oj) . (4) 

The operators with n < 4 account for the dependence of the NN force on the total 
spin and isospin of the interacting pair, while produces the non central interaction 
responsible for the nonvanishing electric quadrupole moment of the deuteron. Inclusion 
of the the operators with n < 6 is needed to obtain a reasonable description of the 
two-nucleon bound state. The operators corresponding to n = 7, . . . , 14 are associated 
with the non static components of the NN interaction, while those corresponding to 
p = 15, . . . , 18 account for charge symmetry violations. The minimal set required to 
describe NN scattering in S and P states (i.e. states of relative angular momentum 

7 8 

L = 0, 1) consists of the operators corresponding to n < 8, with = L ■ S® [1, (t,- • T/)]. 

The ellipses in Eq.© refer to the presence of interactions involving more than two 
nucleons. It is long known that the inclusion of a three-nucleon potential Vijk is needed 
to reproduce the observed binding energies of the three-nucleon system, as well as 
the equilibrium density of symmetric nuclear matter [5]. Note that, while at nuclear 
density CLk>j>i^ijk) <^ (Lj>i v ij)> m tne density regime relevant to the description of 
the neutron- star core many body forces are expected to provide a sizable contribution to 
the energy. 



Correlated Basis Functions 



Due to the presence of a strongly repulsive core, the matrix elements of the NN inter- 
action between eigenstates of the non interacting system (Fermi gas (FG) states in nu- 
clear matter) are large. As a consequence, standard perturbation theory is not applicable. 
This difficulty can be overcome carrying out a resummation of the perturbative expan- 
sion, leading to the replacement of the bare interaction with the well behaved G-matrix 
(for a recent review see, e.g., Ref. [6]). An alternative approach, originally developed to 
study strongly correlated quantum liquids [Q], is based on the replacement of the FG ba- 
sis with a basis of correlated states, embodying the nonperturbative effects arising from 
the short range behavior of the NN force. 

The correlated states of nuclear matter are obtained from FG states through the 
transformation 

\n)=F\n FG ), (5) 

where the operator F, carrying the correlation structure induced by the NN interaction, 
is written in the form 

F= ^n^' (6) 

ij 

y being the symmetrization operator. The two-body correlation functions fy, whose 
operatorial structure reflects the complexity of the NN potential, read 

/y=E/"(ry)0&, (7) 
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with the O'-j given by Eq.©. 

The radial functions f n (rij) are obtained from functional minimization of the expec- 
tation value of the nuclear hamiltonian (QQ) in the correlated ground state 



The effective interaction 

The effective interaction v e ff is defined by the relation [Q]] 

(H) = ^^ = (0 FG \K + v eif \0 FG ), (8) 

where H is the full nuclear hamiltonian of Eq.© and K is the kinetic energy operator. 

In order to include interactions involving more than two nucleons, we have followed 
the approach originally proposed in Ref. [9], in which the main effect of three- and 
many-body forces is taken into account through a density dependent modification of the 
NN potential v,; ,■ at intermediate range. Moreover, in view of the weak model dependence 
of our results [2], the full vis potential has been replaced by its reduced form v 8 lllOll . in 
which only the components with n < 8 are retained. 

The cluster expansion technique LI 1] allows one to rewrite (H) in the form 

(H)=K FG +Y,(*E) n , (9) 

n>2 



where KfQ is the Fermi gas kinetic energy and (AE)„ is the contribution to the energy 
arising from n-nucleon clusters. 

To obtain v e ff from Eq.©, (H) is evaluated at the two-body level of the cluster 
expansion JH]. The resulting effective interaction reads 
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For any given density, the radial functions / n (r ? j) of Eq.© are solutions of a 
set of Euler-Lagrange equations satisfying the boundary conditions f^irtj > d) = 1, 
f n ( r ij > d) = 0, for n = 2, 3 and 4, and f n {nj > d t ) = 0, for i = 5, 6 (see, e.g., Ref. 
[9]). Note that, as the non static terms in Eq. (fTUl) yield a negligibly small contribution to 
the ground state energy yj], they have been neglected altogether. 

The effective interaction of Eq. (fT0l) was tested by computing the energy per particle 
of symmetric nuclear matter and pure neutron matter in first order perturbation theory 
using the FG basis. In Fig. [T]our results are compared to those of Refs. [8] and [12]. The 
calculations of Ref. (SO (solid lines) have been carried out using a variational approach 
based on the FHNC-SOC formalism, with a hamiltonian including the Argonne vis NN 
potential and the Urbana IX three-body potential [5]. The results of Ref. (dashed 
line of the lower panel) have been obtained using the Vg and the same three-body 
potential within the framework of the Auxiliary Field Diffusion Monte Carlo (AFDMC) 
approach. 
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FIGURE 1. Energy per particle of symmetric nuclear matter (upper panel) and pure neutron matter 
(lower panel). The diamonds represent the results obtained using the effective interaction discussed in the 
text in first order perturbation theory with the FG basis, whereas the solid lines correspond to the results 
of Akmal, Pandharipande and Ravenhall [8]. The dashed line of the lower panel represents the results of 
the AFDMC approach or Ref. ifll . 



The results of Fig. Q] show that the effective interaction provides a fairly reasonable 
description of the EOS over a broad density range. Note that empirical equilibrium prop- 



erties of symmetric nuclear matter are accounted for without including the somewhat ad 
hoc density dependent correction of Ref. [8]. It should also be emphasized that our ap- 
proach does not involve adjustable parameters. The correlation ranges d and d t have been 
taken from Ref. (SD, while the parameters entering the definition of the three-nucleon in- 
teraction (TNI) have been determined by the authors of Ref. fl9J] through a fit of nuclear 
matter equilibrium properties. 



Transport properties of interacting Fermi systems 

The theoretical description of transport properties of normal Fermi liquids is based on 
Landau theory f| 1 3T | - Working within this framework and including the leading term in 
the low-temperature expansion, Abrikosov and Khalatnikov (AK) [14] obtained approx- 
imate expressions for the shear viscosity and the thermal conductivity. Let us consider 
viscosity, as an example. The AK result reads 

where p is the density, vp is the Fermi velocity and m* and x denote the quasiparticle 
effective mass and lifetime, respectively. The latter can be written in terms of the angle- 
averaged scattering probability, (W), according to 



m* 3 {W) ' 



where T is the temperature and 



{w) = j2^^{ej2)- (13) 

Note that the scattering process involves quasiparticles on the Fermi surface. As a 
consequence, for any given density p, W depends only on the angular variables 6 and 
0, the magnitude of all quasiparticle momenta being equal to the Fermi momentum 
Pf = (37T 2 p) 1//3 . Finally, the quantity A.^ appearing in Eq. tfTTb is defined as 

(W[l-3sin 4 (6/2)sin 2 0]) 

77 = W) ■ ( } 



The exact solution of the equation derived in Ref. II 1411 . obtained by Brooker and Sykes 
UM, reads 

1-/L~ 4k + 3 

1 = r\AK —± £ 



4 ^(*+l)(2*+l)[(*+l)(2*+l)-A^] 



the size of the correction with respect to the result of Eq.(fm) being 0.750 <(f\f t\ak) < 
0.925. 



Eqs. (fTTI) -([T5l) show that the key element in the determination of the viscosity is the 
in-medium NN scattering cross section. In Ref. Rial , the relation between NN scattering 
in vacuum and in nuclear matter has been analyzed under the assumption that the nuclear 
medium mainly affects the flux of incoming particles and the phase space available to 
the final state particles, while leaving the transition probability unchanged. Within this 
picture W(6, 0) can be extracted from the NN scattering cross section measured in free 
space, (do / dQ.) vac , according to 
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do 
dQ. 



(15) 



where m* is the nucleon effective mass and 6 and are related to the kinematical 
variables in the center of mass frame through E cm = p\ ( 1 — cos ) / (2m) , cm = (j) (vft\ . 

The above procedure has been followed in Ref. II 1811 . whose authors have used the 
available tables of vacuum cross sections obtained from partial wave analysis ||19I1 . In 
order to compare with the results of Ref. [18], and gauge the model dependence of our 
results, we have carried out a calculation of the viscosity using Eqs. ((TT|) -(IT5l) and the free 
space neutron-neutron cross section obtained from the Argonne vis and v 8 potentials. 
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FIGURE 2. Differential neutron-neutron scattering cross section at E cm = 100 MeV, as a function of 
the scatterin g an gle in the center of mass frame. Solid line: cross section in vacuum, calculated with the 
v' g potential IllOtl . Dot-dash line: medium modified cross section obtained from the effective interaction 
described in the text at p = 0.08 fm~ 3 . Dashed line: same as the dot-dash line, but for p =0.16 fm~ 3 . 

In the upper left panel of Fig. [3] we show the quantity r]T 2 as a function of density. 
Our results are represented by the solid line, while the dot-dash line corresponds to the 
results obtained from Eqs.(43) and (46) of Ref. [18] using the same effective masses, 
computed from the effective interaction discussed above. The differences between the 
two curves are likely to be ascribed to the to the extrapolation needed to determine the 
cross sections at small angles within the approach of Ref. f| 1 8T | - 

To improve upon the approximation of Eq.dTSl) and include the effects of medium 
modifications of the NN scattering amplitude, we have replaced the bare NN potential 
with the CBF effective interaction. 

Knowing the effective interaction, the in-medium scattering probability can be readily 
obtained from Fermi's golden rule. The corresponding cross section at momentum 



transfer q reads 



(16) 



v e ff being the Fourier transform of the effective potential. The effective mass can also be 
extracted from the quasiparticle energies computed in Hartree-Fock approximation. For 
symmetric nuclear matter at equilibrium, we find m*(pf)/m = 0.65, in close agreement 
with the lowest order CBF result of Ref. [200 . 

In Fig. |2]the in-medium neutron-neutron cross section at E cm = 100 MeV obtained 
from the effective potential, with p = po and po/2 (po =0.16 fm~ 3 is the equilibrium 
density of nuclear matter), is compared to the corresponding free space result. As 
expected, screening of the bare interaction leads to an appreciable suppression of the 
scattering cross section. 
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FIGURE 3. Left Panel: neutron matter 7]T 2 as a function of density. (A). Solid line: results obtained 
from Eqs.dTTb-dTSb using the Argonne Vig potential and m* computed from the effective interaction 
described in the text. Dot-dash line: results obtained from Eqs.(43) and (46) of Ref. [18] using the same 
m*. Dashed line: same as the solid line, but with the Argonne vis replaced by its reduced form v' 8 . (B). 
Solid line: results obtained using the effective interaction described in the text. Dashed line: r\T 2 obtained 
from the free space cross section corresponding to the v' 8 potential. Right panel: neutron matter kT as a 
function of density. Solid line: calculation carried out using the approach of Ref. [2]; dashed line: results 
of Baiko, Hansel & Yakovlev 12211 . 

Replacing the cross section in vacuum with the one defined in Eq.(fT6l), the medium 
modified scattering probability can be obtained from Eq.(fT5l). The resulting W(0, 0) can 
then be used to calculate rjT 2 from Eqs. (fTT]) -(fT5T). Similar expression are found for the 
thermal conductivity. 

The effect of using the medium modified cross section is illustrated in the lower left 
panel of Fig. [3l Comparison between the solid and dashed lines shows that inclusion of 
medium modifications leads to a large increase of the viscosity, ranging between a factor 
~ 2.5 at half nuclear matter density to a factor ~ 10 at p = 2po- Such an increase is likely 



to produce appreciable effects on the damping of neutron-star oscillations associated 
with emission of gravitational waves 112 ill . 

The right panel of figure [3] shows that medium modification of NN scattering also 
produce a large effect on the thermal conductivity, K. The comparison between our 
results and those of Ref. [22], obtained from the free space NN cross sections, shows 
a difference in kT exceeding one order of magnitude at p = 2pQ. The impact of this 
result on neutron star cooling is currently being investigated |23|]. 



SUMMARY 

Working within the formalism based on correlated functions and cluster expansion 
techniques, we have derived an effective interaction from a realistic NN potential model. 
Our work improves upon the CBF effective interaction of Ref. yj] in that it includes the 
effects of many-nucleon forces, which become sizable, indeed dominant, in the high 
density region relevant to the studies of neutron star properties. The energy per nucleon 
of both symmetric nuclear matter and pure neutron matter, obtained from our effective 
interaction model, turns out to be in fairly good agreement with the results of highly 
refined many body calculations, based on similar dynamical models. The emerging 
picture suggests that our approach captures the relevant physics, allowing for a unified 
description of a number of different properties of neutron star matter based on standard 
perturbation theory in the FG basis. 
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